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VIEWED AS AN APPLICATION OF THE POYNTING THEOREM 
by John Barranger 
Lewis Research Center 

SUMMARY 

An energy -loss equation is developed based on the Poynting theorem. The resultant 
integrals provide a unified viewpoint that is shown to be the same as the classical eddy- 
current and hysteresis energy formulations. An expression is derived for the eddy- 
current loss of a thin transformer lamination by assuming that the permeability is con- 
stant, and further, that the magnetic material is isotropic and homogeneous. The hys- 
teresis loss is found by assuming that the introduction of the hysteresis loop does not 
alter the fields. The results indicate that, for a given total flux and at frequencies ap- 
proaching total shielding, both losses are proportional to the 3/2 power of frequency and 
the square root of the conductivity and inversely proportional to the square root of 
permeability. 


INTRODUCTION 

As part of a study of the distribution and use of power on spacecraft, the Lewis 
Research Center has been evaluating power system components in unusual configurations 
and at temperatures and frequencies outside the ordinary ranges. Early in this study it 
became evident that neither classical nor modern approaches to magnetism could com- 
pletely predict transformer properties and their variations with frequency. 

The classical approach to the study of losses in ferromagnetic materials is the con- 
ventional engineering analysis which assumes the magnetic material to be isotropic and 
homogeneous and defines the losses in terms of the hysteresis loop and classical electro- 
dynamics. The modern theory of magnetism utilizes crystal energies and magnetic do- 
main configurations to provide an understanding of losses (refs. 1 to 3). The classically 
derived losses are characterized by measurements that are relatively easy to make while 
the modern theory relies on crystal anisotropy energies and domain wall velocities that 
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are very difficult to evaluate and measure in polycrystalline materials. 

Traditionally, transformer core losses have been divided into two groups, hysteresis 
and eddy-current loss. Hysteresis loss is thought to be associated with the energy re- 
quired to rotate or move the walls of the magnetic domains over a full cycle. The hys- 
teresis loss is proportional to the area of the familiar hysteresis loop. The eddy-current 
loss is due to the currents generated within a real conductor subjected to a varying mag- 
netic field. 

The purpose of this report is (1) to provide a better understanding of classical theory 
for the readers already familiar with the basic engineering approach and (2) to derive ex- 
pressions useful in predicting the dependence of core losses on frequency and easily 
measured properties. Hysteresis and eddy-current losses are usually treated separately 
with no common underlying basis. A more unified viewpoint will be provided starting 
with Maxwell’s equations and will show that considerations of the Poynting theorem yield 
the classical core loss formulations. The separation of core losses into their two tra- 
ditional components results as a natural consequence of the analysis. 

A large part of the study is devoted to finding expressions for the hysteresis and 
eddy-current losses for thin transformer laminations. The magnetic material will be 
assumed to be isotropic and homogeneous and to have a permeability that is constant. 

The permeability is obtained from the mean magnetization curve by taking the ratio of the 
flux density to the corresponding magnetic field intensity. The formulas will then be ex- 
amined to ascertain the behavior of the losses with frequency. To simplify the resultant 
expressions, low- and high-frequency approximations will be derived. The results will 
also be used to predict the important parameters of use in high-frequency transformer 
design. This should be of interest to designers who are seeking methods of using lami- 
nated magnetic material for high-frequency power applications. It is possible that im- 
provements in weight and power efficiencies may be gained through the use of high- 
frequency power systems. 


SYMBOLS 

a lamination thickness 
B flux density 
b length of lamination 
E potential 
f frequency 
H field intensity 
J current density 
2 


Z height of lamination 
P power loss 

5 shape factor 

t time variable 

V volume 
W energy 

6 skin depth 

6 dummy variable 
M permeability 
p free charge density 
a conductivity 
$ total flux 


POYNTING THEOREM AND CORE LOSSES 

Energy transfer from a source to a magnetic material is accomplished through elec- 
tromagnetic induction, the amount of energy depending on the magnitude, distribution, 
and phases of the electric and magnetic fields. Maxwell's equations, written in terms of 
the total fields and currents within the magnetic material, describe the electromagnetic 
behavior of the material. The following conditions are appropriate to these equations 
when applied to good conductors: 

(1) The free-charge term is zero, p = 0. __ 

(2) Conduction current density is given by Ohm's law, J(t) = aE(t). 

(3) Displacement current is negligible in comparison with conduction current. 

(4) The material is isotropic. 

The two curl equations become 


V X E(t) = - 

at 


(la) 


V X H(t) = J(t) 

If the divergence of E(t) x H(t) (the Poynting vector) is considered, 


(lb) 
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V • [i(t) X H(t)] = H(t) • [v X E(t)] - E(t) • [v X H(t)] = - H(t) . ML - E(t) • J(t) 

3t 


(2) 


Integrating over the enclosed volume produces 


f V • [E(t) x H(t)]dV = - J |^H(t) . ML + E(t) • J(t)JdV 


(3) 


and using the divergence theorem results in 


- J [E(t) x H(t)] • dS = J ^H(t) • MD + E(t) • J(t)jdV 


(4) 


The total of all the terms to the right of the equal sign are identified as representing the 
energy stored and dissipated per unit time, the left side of equation (4) represents the 
energy flow into the volume per unit time. The energy flow into the volume, then, be- 
tween two points of time t^ and tg is 


W = 



3B(t) 

at 


+ E(t) • J(t) 


dVdt 


(5) 


The time rate of change of flux density for a stationary volume is 


dB(t) _ 8B(t) 


dt 


at 


and 


W 


= f 2 / [n(t) . MD + E(t) • j(t) 
J. J-tr L dt 


'tj "V 


dV dt 


( 6 ) 

(7) 


Since the material was assumed isotropic, B(t) and H(t), E(t) and J(t) are in the same 
direction and the energy expression becomes 
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where Bj and Bg are the flux densities corresponding to ^ and t 2 , respectively. 
The second term represents the energy loss due to conduction currents or the eddy- 
current loss. Over a complete cycle, the first term represents the energy loss, which 
is called hysteresis loss. Thus, according to Maxwell’s equations the total loss or core 
loss is 


P c = P h + P e < 9 > 

where P^ is the hysteresis loss and P 0 is the eddy-current loss. 

Consider now the nonlinearities commonly encountered in ferromagnetic substances. 
The two significant nonlinearities are 

(1) The hysteresis loop 

(2) The mean magnetization curve 

To simplify the analysis, the mean magnetization curve will be assumed to be a straight 
line, that is, the permeability is a constant. The material will also be assumed to be 
homogeneous as well as isotropic. 


CALCULATION OF CURRENT AND FIELD 

Before any calculation of core loss can be made, the current and field distribution 
within the volume of the material must be known. This is usually accomplished by solving 
a differential equation in the quantity of interest. By returning once more to Maxwell’s 
equations, the curl of the magnetizing force H(t) may be written 

V x H(t) = <rE(t) (10) 

and taking the curl of both sides yields 

-V 2 H(t) + V [v • H(t)] = -a (11) 

3t 

Since it is desirable to solve the previous expression in terms of one variable, such as 
B(t), the classic problem of relating B(t) and H(t) for ferromagnetic substances must be 
faced. Even though the hysteresis loop shows that this is a nonlinear functional relation, 
it will be assumed that the introduction of the loop alters only slightly the distribution 
and waveshape of the fields . Oriented silicon iron approximates this condition at low to 
moderate flux density levels. The functional relation then becomes 
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B(t) = /xH(t) 


( 12 ) 


where the permeability p. is the constant slope of the mean 
magnetization curve. Substituting, the resultant equation 
becomes 


V 2 B(t) = CTjU 

at 


(13) 


Consider the lamination of a transformer as a thin sheet 
of thickness a where the height of the lamination l is much 
greater than the thickness. Let the field be applied to the sur- 
face of the lamination in such a way that the induced eddy cur- 
rents are as indicated in figure 1. It is necessary to know how these currents are dis- 
tributed before an expression for power loss can be found. 

The equation derived previously for the variation of B(t) in a conductor is in the form 
of a standard differential equation similar to the diffusion or heat equation. By a similar 
derivation, the equation that gives the relation between space and time derivatives of 
current density at any point in a conductor is 


V 2 J(t) = <TM 

at 


(14) 


Letting J(t) be sinusoidal and introducing the complex exponential form 

J(t) = Re(je jwt ) 

results in the previous expression becoming 


(15) 


V J = JcdctmJ 


(16) 


Study of figure 1 will show that, since no net current can flow out of a conductor, 
the following symmetry and boundary conditions exist for thin sheets with an axially 
applied field: 

(1) J x = 0 


(2) J y = 0 
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(3) 


9J„ _ 9% 


= 0 


(4) 


3 j„ 


9z 9z 2 


= 0 


(5) J, 


(6) <L 


x=(a/2) J ° 


x=-(a/2) 


= -J. 


By using the first four conditions, the current distribution equation becomes simply 


d 2 J 


- = jcoajuJ z = T 2 J z 


(17) 


dx 


where T is a constant and is equal to (taking the root with positive sign) 


T = (1 + j)-^7rf j Ha 


(18) 


It is convenient to write the solution of the differential equation in terms of hyperbolic 
functions 


J = A sinh Tx + C cosh Tx 
z 


(19) 


To evaluate arbitrary constants A and C, boundary conditions (5) and (6) (above) are 
utilized 


J = A sinh T - + C cosh T - 
2 2 


- J = -A sinh T - + C cosh T — 
° 2 2 


(20a) 

(20b) 


Solving these results in 


C = 0 


(21a) 
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(21b) 


A = 



sinh T - 
2 


and 


J = J 
z o 


sinh Tx 
sinh T - 


(21c) 


Defining the quantity 6 as \jyfrd\Io, the current density distribution may be written as 

Hjf) 


j = j 

Z 0 




where 6 is usually called the skin depth or depth of penetration. In terms of a semi- 
infinite solid, it is the depth at which current density has decreased to 1/e (about 36. 9 
percent) of its value at the surface. There are a number of important factors associated 
with the skin depth; they are the following: 

(1) The depth of penetration decreases as the conductivity, the permeability and the 
frequency increases. Thus, materials of high permeability have a smaller depth of 
penetration, or conversely, skin effect takes place at a lower frequency. 

(2) The current does not fail to penetrate below the depth 6; this is merely the point 
at which current densities and fields have decreased to 1/e of their value at the surface. 

(3) The skin depth may be considered a constant for a given material at a frequency 
f and is a useful parameter in the analysis of thin sheets. 

The current distribution at low frequencies may be obtained by making the ratio 
a/26 much less than one, that is, the skin depth is much greater than the thickness. 
Taking the absolute value of the current density gives 


Jj = |J 


sinh — cos — + j cosh - sin — 


6 


6 

cl 3L 3. a. 

sinh — cos — + j cosh — sin — 

26 26 26 26 



= J. 


2 x 2 x 
sinh^ — + sin z — 


. ,2 a .2a 
sinh — + sin — 


(23) 


26 


26 


and if a/26 « 1 and x/6 « 1, 
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I 




(24) 


This distribution is linear with x. There remains only Jj Q | to be evaluated. 

Using Faraday's line integral and integrating around the path ABCD in figure 1 yield 



(25) 


where «f>(t) is the total flux within the line ABCD. Multiplying both sides by conductivity 
a and using the relation J(t) Q = ctE ( t) Q make the following obtainable: 


/j(t) 

J ° dt 


(26) 


If the sheet has height l and the height is much greater than the thickness, the current 
density J Q in complex notation form is 


If the absolute value is taken 


then 


■ _ jh)cr$ 
0 21 



i oo |$ 
21 



. . 2 x , .2 x 

sinh — + sin — 


1/2 


..2a .2a 

, sinh — + sin — 


26 


26. 


(27) 


(28) 


(29) 


The previous expression represents the current density distribution, which will be 
used to find the eddy-current loss. In the analysis of hysteresis loss, the flux density 
distribution must also be known. This may be found directly from the current density 
expression. Consider the curl equation 


9 


I 
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V x E(t) = - d-ft) 

at 


(30) 


Using the complex exponential form 


V x J = - jo>aB 


( 31 ) 


and noting that J is the only component of J result in 
z 


dJ z 

= - jcocrB 

0X y 


( 32 ) 


or 


i 9J z 

B y = ' — 

J u)<y dx 


( 33 ) 


If the formulas for J z and J Q found previously are now substituted, 


V 


_ j 9_ /j sinh Tx \ _ #T / cosh Tx 

coa 9x \ ° sinh T a I 21 ( sinh T a 


( 34 ) 


The absolute value of the flux density becomes 


/cos 2 i+ Sinh 2 5 \‘ /2 

= Vgj$! I 6 6 1 


y' 2 i 6 


i . 2 a . ,2 a 
(sin — + sinh — 

\ 26 26> 


( 35 ) 


EDDY-CURRENT LOSS 


The first section showed that the eddy-current loss may be expressed as 


ft 

V •'tj 


2 dt dV 


( 36 ) 
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It can be shown that the average power loss in complex notation is 

|2 


P e = 


I 


2a 


dV 


( 37 ) 


If the element of volume is bZ dx (fig. 2), then the total power may be expressed as 


P e = 2 


/ 


a/ 2 


z ' bZ dx = — 


2a 


\ - f (sinh 2 — + sin 2 — \ dx 

' sinh 2 — + sin 2 — *1) V 6 


25 


25 


co 2 |#| 2 ab5 

81 


sinh — - sin — 

5 5 

, a a 

. cosh — cos — 

6 


(38) 


The previous expression represents the power loss in a transformer lamination due 

to eddy currents. Latour (ref. 4) considered the same problem from the standpoint of a 

constant lag of phase angle between the magnetic induction B(t) and magnetizing field 

H(t). His loss formula is the same as the one derived (eq. 38). For a given lamination 

o 

material, thickness, and frequency, the loss varies as <£ . Simpler expressions may 

be obtained by examing the behavior of the equation as a func- 
tion of frequency. At low frequencies or when a/5 « 1, the 
power loss is 



e(L) 


o) 2 4> 2 aba 
max 

24 1 


(39) 


where is the maximum value of total flux. 

niaX 


If B is 

assumed to be uniform at power frequencies, the average 
power per unit volume yields the familiar result 


5 e(L) 


• 2 ‘ 2 4a/ V 


(40) 


Figure Z - Lamination of transformer 
showing element of volume. 


At frequencies where skin effect predominates, that is, 
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when a/ 6 » 1, the power loss is 


P e(H) - 


crb6 

max 


81 


3/2 f 3/2^2 l/2 b 
max 

2 1 M 1 ^ 2 


(41) 


and the average power per unit volume is 


p e(H) - 


„3/ V /2 # 2 1/2 


max 


2l 2 aM 1/2 


(42) 


Thus, if the sinusoidal field is regulated so that there is constant $ max , the eddy- 
current loss varies as the square of the frequency at power frequencies and as the 3/2 
power of frequency at frequencies approaching total eddy-current shielding. 


HYSTERESIS LOSS 


It was stated previously that the integral 



H dB dV 


B 



represents the hysteresis energy loss when it 
is taken over a full cycle, that is, 


W(per cycle) = 



dV 


(43) 


The bracket represents the area of the hyster- 
esis loop for the element dV under considera- 
tion. If the area is expressed as a function of 
the variable V, then the formula may be inte- 
grated over the entire volume to yield the total 
hysteresis loss per cycle. The H and B 
within the bracket are obtained from the static 
hysteresis loop of the material. The end point 
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of the hysteresis loop corresponds to the maximum value of B. Thus, the solution of the 
differential equation in B may also be used to describe this point. The area of the 
hysteresis loop will be found first for the surface of the sheet, and then the results will 
be extended to provide an expression for the area at any point. 

The hysteresis loop for the surface is given in figure 3. To find the area of this 
loop, a dummy variable 6 defined by the following equation must be introduced: 

B = -B cos 6 (44) 

where B m is the maximum value of B. If 6 is limited to a half period, the magnetiz- 
ing force H, may be described by the half -range Fourier expansion 

H = a n H + a,H, cos 9 + a 0 H, cos 26 + • • • + a„H, cos n 6 (45) 

Umlm zm nm v 

where H m is the maximum value of H. Examination of figure 3 shows that, as the flux 
density traverses the loop from -B to +B m and back to -B , the area over the 
complete cycle is twice the area between -B to +B m . Thus, 



H dB = 2 



H dB 


(46) 


and if it is noted that 6 = 0 and 6 = u correspond to -B and +B m respectively, the 
area becomes 


= 2 f 

Jr \ 


HB 
m m 


sin 6 


( a 0 


+ a« cos 6 + 


cos 


26 + 


a n cos nfljde 


= 2H m B m 






for n even 



for n even 


(47) 


Let the summation within the parentheses be called the shape factor S. The value of the 
shape factor may be determined by solving the parametric equations for B and H for 
the hysteresis loop being investigated. The number of terms in the summation depends 
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on the number of data points and the closeness of fit desired. Since the permeability was 
assumed constant, B m and H m are related by 


B„ = /ill 
m m 


(48) 


and the area of the hysteresis loop is 2B^S j / u. 

The hysteresis energy integral can now be solved 


W(per cycle) = / [/(HdB)]d V=y 


2B S 
m 


dV 


(49) 


According to figure 2 (p. 11), 


W(per cycle = 4 lb 


/»a/2 

J n 


B 2 S 
m 


dx 


(50) 


and since S and ju are assumed constant, they will not depend on the coordinate x. 
Thus, 


W(per cycle) = 


4lbS 


f a/2 
*'0 


m 


(51) 


It is now necessary to find B in terms of x to integrate the expression. It was found 
previously that 


+ sinh 2 -\^ 2 


bJ = 



(52) 


and by noting that 


B = B 
y m 


( 53 ) 



may now be solved 
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I 


B 2 -2|± 
m 


„„„2 x , • ,2 x 

cos — + sinn — 


(2Z6) 2 Isin 2 — + sinh 2 — 


25 


26 , 


(54) 


Substituting into the energy equation yields 


W(per cycle) = 


2 IbS 
M 


m ) 2 — _! r /2 (, 

\l 6 / ..2a ■ 2 a. J n V 

' ' sinh — + sin — •'O 


cos 2 — + sinh 2 — )dx 


26 


26 


r sinh — + sin — 

bsi$r 

\cosh — - cos — , 


(55) 


The hysteresis loss is f times the energy lost per cycle. Thus, 


P h = fw = 


fbS 


F sinh — + sin — 


, a a 

icosh — - cos — 


(56) 


The previous expression represents the power loss in a transformer lamination due 
to hysteresis effects. Latour (ref. 4) found the same formula except that this formula 
for hysteresis loss has the added advantage of being defined in terms of data taken from 
the static hysteresis loop. At power frequencies or when a/6 « 1, the power loss is 


P h(L) - 


2fbS#: 


max 


lil a 


(57) 


If B is assumed uniform at power frequencies, the average power per unit volume 
yields the result 


p h(L) - 


2fSB 


max 




(58) 
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which again is frequency times the area of the hysteresis loop. 

At frequencies where skin effect predominates (i.e. , when a/6 » 1) the power loss 
is 


fbS* 2 a 1 / 2 

p h(H) - 


and the average power per .unit volume is 


*>h(H) = 


, 1 / 2 , 3 / 28*2 1/2 

m ax 

_ p 1 / 2 1 2 a 


(59) 


(60) 


Thus, if the sinusoidal field is regulated so that there is constant 
as f 2 / 2 . 


<$. the loss varies 

max’ 


TOTAL CORE LOSS 


At high frequencies, the hysteresis and eddy-current losses can be combined to 
form the total core losses 


p c(H) " p e(H) + p h(H) 


, tSB max a 
2 ju6 


( 61 ) 


where the average flux density B max is defined as 


B 


$ 


max 


max 


l a 


(62) 


After substitution, the previous equation yields 


Jr V2,3/2g 2 ,1/2 


■0(H) 


max 


P 


1/2 


(H 


(63) 


For a given B max and frequency, then, the core loss at high frequencies is proportional 
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to the square root of conductivity and inversely proportional to the square root of the 
permeability. Clearly, as the frequency is increased, the eddy-current loss varies 
from the square to the 3/2 power of frequency, while the hysteresis loss varies from the 
first to the 3/2 power of frequency. 


CONCLUDING REMARKS 

It has been demonstrated that the use of Maxwell's equations and the Poynting 
theorem yield the classical eddy-current and hysteresis loss formulations. The resultant 
integrals were evaluated to find expressions for each loss for a thin transformer lami- 
nation by assuming that the permeability was constant, and further, that the magnetic 
material was isotropic and homogeneous. The expressions were then simplified to yield 
low-frequency formulas. The hysteresis loss was found to be proportional to the area of 
the hysteresis loop, while the eddy-current loss became the conventional engineering ex- 
pression. 

By making a high-frequency approximation, it has been found that for a given average 
B max anc * frequency, the core loss at high frequencies is proportional to the conductivity 
and inversely proportional to the permeability. Thus, to minimize losses, high- 
frequency transformer laminations should be made of a material that has high resistivity 
and high permeability. 

Examination of the core loss expressions also showed that for a given lamination 
material, thickness, and $ , as the frequency is increased, the eddy-current loss 

varies from the square to the 3/2 power of frequency while the hysteresis loss varies 
from the first to the 3/2 power of frequency. This means that the eddy-current loss con- 
tribution to the total loss is less than would be expected from the power frequency ex- 
pression and that the hysteresis-loss contribution to the total loss is greater than would 
be expected from the power -frequency considerations alone. Thus, hysteresis loss 
plays an important role in transformer core loss at high as well as low frequencies. 

Since the expressions for the thin transformer lamination assumed that the perme- 
ability was constant, application of the loss formulas without regard to variations of 
permeability with the applied field and saturation effects will result in large errors. It 
is felt that by maintaining the same basic approach, an approximation scheme could be 
devised that would take the true permeability into account (refs. 5 to 7). The assumption 
that the material is isotropic and homogeneous would still have to be made. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 1, 1965. 
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